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ABSTRACT This paper describes a new research effort to 
extend the application of Global Positioning System (GPS) 
signal reflections, received by airborne instruments, to 
cryospheric remote sensing. Our experimental results 
indicate that reflected GPS signals have potential to provide 
information on the presence and condition of sea and fresh- 
water ice as well as the freeze/thaw state of frozen ground. 
In this paper we show results from aircraft experiments over 
the ice pack near Barrow, Alaska indicating correlation 
between forward- scattered GPS returns and RADARSAT 
backscattered measurements. 

INTRODUCTION 

In 1996, Drs. Katzberg and Garrison of NASA Langley 
Research Center (LaRC) developed the idea of using 
reflected GPS signals for remote sensing applications, and 
published a number of papers describing the theory and 
mechanisms for this technique [e.g., 1-3]. Since then, 
research has advanced our understanding of reflected GPS 
signals and experimentally applied those techniques to ocean 
remote sensing and mapping [4-10]. The current investigation 
extends this work into cryospheric applications of GPS 
reflections based upon a proposal made by Drs. Katzberg and 
Garrison, to collect GPS signals reflected by ice surfaces. 
Results are shown from the first GPS measurements from sea 
ice in the Beaufort Sea in April 1998. 

MEASUREMENT METHODS 

The use of GPS in a bistatic radar configuration to measure 
surface properties relies upon our ability to extract 
information from the reflected signal. For standard GPS 
navigation purposes, the GPS receiver measures the signal 
delay from the satellite (the pseudorange measurement) and 
the rate of change of the range (the Doppler measurement). 
For remote sensing applications, the primary measurement is 
the received power from the reflected signal for a number of 
delays and Doppler values. This measurement and its 
sensitivity to the surface conditions are discussed in the 
following paragraphs. 

The GPS satellites transmit two right-hand circularly 
polarized (RHCP) L-band signals at 1.57542 GHz and 1.2276 
GHz. These carrier signals are modulated by unique 
pseudorandom noise (PRN) codes with an autocorrelation 
function similar to the ideal autocorrelation function. The 
autocorrelation power has a triangular shape, which is 


suppressed at delays of more than ± 1 chip (300 m for the 
civilian (C/A) code). To acquire and track a signal, a 
conventional GPS receiver generates a local replica of the 
code for the particular satellite signal to be tracked, 
compensating for the expected Doppler shift and computing 
the cross-correlation between the replica and the incoming 
codes. 

When a GPS signal encounters an ideal, smooth reflecting 
surface, specular reflection occurs at a single point. The 
reflected signal code structure remains, but the polarization of 
the wave is reversed to left-hand circular polarization 
(LHCP), and the signal power is decreased. If, however, the 
surface is rough relative to the GPS wavelength of 19 cm, 
reflections are produced by multiple facets on the surface. 
Around the ideal specular point, this creates a so-called 
glistening zone within which there is a distribution of varying 
ranges and Doppler shifts, as shown in Fig. 1. To measure 
GPS signals reflected from land and sea surfaces, Garrison 
and Katzberg [1] modified a typical receiver to measure 
correlation power at these offset values of delay and Doppler 
using a nadir-pointing LHCP antenna. 

In Fig. 2, we illustrate the type of swath coverage that GPS 
provides. Because of the multiplicity of GPS transmitters the 
receiver can observe a number of simultaneous footprints on 
the surface. This illustrates that while the system does not 
provide imaging in a standard sense (like a cross-scanning 
system, for example), it does provide more spatial coverage 
than from a sensor that is only able to view one ground point 
at a time. 



Figure 1. Illustration of glistening zone (adapted from 
Komjathy etal., 1999). 
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Figure 2. GPS swath coverage. 

Results and observations suggest that the reflected GPS 
signal is very sensitive to the presence of sea ice, and 
furthermore, contains information related to the reflection 
coefficient of different ice conditions. 

The reflection coefficient of a frozen sea surface is 
determined by the effective dielectric constant of ice and, 
under some conditions, by the dielectric constant of the 
underlying water. The latter is important for thin ice 
conditions where a significant portion of the radio wave 
energy may reach the second interface, between the ice and 
water. The effective dielectric constant of ice depends on 
various factors, such as an ice composition, salinity, 
temperature, density, age, origin, morphology, etc. [11]. 
Because first-year ice has needle-like inclusions of brine 
(predominantly oriented along the vertical direction), this 
variety of ice is notably anisotropic. Therefore, the dielectric 
constant is a tensor rather than a scalar constant. This makes 
the problem of modeling and interpreting the GPS signal 
scattering from ice rather challenging. At the same time, the 
sensitivity of the GPS scattered signal to these complex ice 
characteristics makes it a particularly attractive remote 
sensing tool. 

Given the unique potential offered by reflected GPS signal 
sensing (e.g., L-band observations, forward-scattering 
sensing, low-cost and simple devices useable on any type of 
aircraft, etc.) and the geophysical significance of additional 
sea ice and permafrost data, research to investigate and 
exploit GPS returns from ice and frozen ground is warranted. 
However much remains to be done to determine the 
consistency of the signal return from ice surfaces; to 
investigate the information content of the signals; and to 
develop simplified retrieval algorithms to extract the 
cryospheric conditions. 

EXPERIMENTS AND DATA ANALYSIS 

Most of the GPS reflection experiments to date have been 
conducted with the Delay Mapping Receiver (DMR) 
designed by Drs. Katzberg and Garrison [1], based on the 
GEC Plessey (now MITEL) GPSBuilder-2 [12]. The DMR 
has two antenna inputs to allow standard tracking of direct 
signals using a zenith-pointed RHCP antenna and correlation 


measurements of reflected signals using a downward-pointed 
LIICP antenna. DMR makes measurements of the correlation 
between the reflected signal and shifted versions of the local 
signal replica. This provides a trace of the reflected signal 
correlation as a function of code delay, thus mapping the 
return from annular regions of the glistening zone. 

The first measurements of reflected GPS data from sea ice 
were collected by Dr. Maslanik in the Beaufort Sea in April 
1998. During the experiment, the GPS system was tested 
with antennas mounted on a boom and held over various ice 
surfaces at different heights ranging from 30 cm to 20 m. 
Each antenna was mounted on a small 20 cm by 20 cm 
aluminum ground plane, with the two plates then attached to 
provide the up-looking and down-looking mount for the 
antennas. A second reflected GPS data set was collected 
north of Barrow, Alaska in April 1999, using the same system 
but in this case mounted on a Cessna 185. Approximately 2 
hours of data were collected in conjunction with a National 
Ice Center (NIC) reconnaissance flight and nearly coincident 
in time with a RADARSAT overpass. 

Our experiment to observe GPS reflections from Arctic sea 
ice showed that for the moderate altitudes of the airborne 
GPS receiver, the reflected signal shape was fairly consistent 
(sharp, narrow waveform throughout the flight) indicating 
that ice surface roughness variations were not significant at 
L-band. On the other hand, the peak power was found to 
change significantly along the flight track. This behavior of 
the signal is a clear indication of the sensitivity to ice 
reflection coefficients. Correlations seen in our comparison 
between RADARSAT backscatter and GPS forward scattered 
data indicate that the GPS signal may indeed provide useful 
information regarding ice conditions, in addition to the basic 
ability to detect the presence of sea ice. 

Fig. 3 shows the peak return power from GPS satellite 
PRN 30, collected during the flight. Reflected power is 
plotted versus RADARSAT backscatter along the flight leg, 
as a function of along track epochs. The two sets of curves 
indicate a positive correlation (see Labels 1 to 6) between the 
forward -scattered GPS peak power and the backscattared 
RADARSAT data for locations of high reflection 
coefficients, as well as some other variations in the GPS data 
that differ from the RADARSAT backscatter. Observation of 
positive rather than negative correlation requires further 
investigation. 


In Fig. 4 we present a photo taken from the aircraft during the 
flight. The picture were taken when the aircraft flew by Label 
5. The correspondence of the reflected peaks with the 
occurrence of large leads in the same area suggests that the 
high signal returns may have originated from patches of new 
or young ice as observed from the aircraft. In Fig. 4, we also 
indicated an illustration of the possible specular reflection 
track crossing a region of recently formed ice resulting in 
high GPS signal returns shown as gray line. Low signal 



returns around this region may have been due to thicker, first 
year ice conditions (dark line). 



Figure 3. Normalized signal power for PRN 30. 



Figure 4. Photo taken from the aircraft in the vicinity of Label 
5 showing black and darker gray areas of newly formed 
leads, new and young ice. 

CONCLUSIONS 

Field experiments, including data acquired from aircraft 
flights over the ice pack near Barrow, Alaska, suggest that the 
reflected GPS signals contain useful information over sea ice. 
Given this new potential application for GPS remote sensing, 
an effort to investigate reflected GPS signals over sea ice is 
discussed in the paper. A combination of modeling 
considerations, in-situ measurements and aircraft 
observations have been presented to quantify the theoretical 
and observed relationships between reflected GPS signals and 
cryospheric conditions. 
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ABSTRACT - Initial estimates of ocean surface wind 
speed and direction based on observations of reflected 
GPS signals are presented and compared to QuikSCAT 
wind fields. The two wind speed estimates are generally 
in agreement to within 2-3 m/s, and under favorable con- 
ditions of well developed seas and stable wands, direction 
estimates agree to within 10 deg. An overview of the GPS 
technique is presented as well as a presentation and dis- 
cussion of these first results. 

INTRODUCTION 

In 1996, NASA researchers Stephen J. Katzberg and 
James L. Garrison identified the potential application of 
the Global Positioning System (GPS) to remote sensing 
of ocean surface conditions [1, 2], Their concept was 
to use GPS in a bistatic radar configuration with the 
GPS satellite transmitting an L-Band spread spectrum 
signal, and the receiver on an aircraft or spacecraft plat- 
form measuring the reflected signal characteristics. Since 
then, Katzberg and Garrison have continued to investi- 
gate the properties of the ocean-reflected GPS signals and 
have developed a specialized Delay-Mapping GPS receiver 
(DMR) to measure the reflected signals [1, 3]. Using the 
DMR and geophysical models to predict the interaction 
of the GPS signals with sea water, researchers have been 
able to estimate wind speeds on the ocean surface with 
an accuracy of 2 m/s [1]. Recently we have developed a 
new approach that extends the GPS-based wind retrieval 
to determination of wind direction as w r ell as wind speed, 
based on measurements from two or more reflected GPS 
signals [3]. To validate the performance of GPS for wind 
vector retrieval we have previously relied on buoy data 
and wind speed estimates from T0PEX\P0SEID0N [4]. 

This paper describes our first comparisons to wind vec- 
tors derived from the SeaWinds scatterometer onboard 
QuikSCAT. Ku-Band scatterometers represent the cur- 
rent state-of-the-art in sea surface wind remote sens- 
ing, providing highly accurate measurements with good 
temporal and spatial resolution on a global scale. The 
QuikSCAT science team provides a variety of wind prod- 
ucts including wind vector maps at 25 km resolution. 
Nominal accuracy for the QuikSCAT wind vectors is 2 
m/s and 20 deg for wind speed and direction, respectively. 
Performance is somewhat worse for wind speeds below 3 
m/s, for grid points within 25 km of land, and for points 
very close to nadir or at large look angles. Algorithms are 
included in the QuikSCAT data reduction that resolve 
wind direction ambiguities and corresponding wind speed 


estimates. QuikSCAT observations provide a snapshot of 
the ocean surface conditions at the time of the satellite 
overflight. The backscatter measurements at Ku-Band 
are sensitive to rapid changes in the surface conditions 
making this a powerful sensor for wind retrieval. Because 
of the substantially different mechanisms for forward scat- 
tered L-Band reflections, the measurement and retrieval 
of wind vectors from GPS may provide a very useful com- 
plement to scatterometer observations in improving global 
wind observations and models. 

GPS-BASED WIND ESTIMATION 

As presented in [1, 2, 3], GPS-based wind retrievals rely 
on the modification of the correlation function for signals 
reflected by rough surfaces. The basic observable is not 
the pseudorange or carrier phase measurement used for 
most GPS applications, but rather a measure of the cor- 
relation of the reflected signal with a delayed replica of 
the PRN code. The distribution of this function over de- 
lays (and Dopplers) provides a mapping of the reflecting 
surface. In particular we use a model developed by Za- 
vorotny and Voronovich (Z-V) [5, 6] for prediction of the 
reflected signal structure. The Z-V model embodies the 
forward-scatter radar equation with the geometric optics 
limit of the Kirchhoff Approximation. The distribution of 
ocean surface slopes is assumed to be Gaussian with vari- 
ances determined by a wave spectra model such as that 
developed by EJfouhaily et al. [7]. These slope statistics 
define the shape and orientation of the glistening zone 
over w r hich GPS signals are reflected toward the receiving 
antenna. The models currently used predict the maxi- 
mum variance in the up/dow r n-wind direction and mini- 
mum variance in the cross-wind direction. This allows for 
the identification of wind direction with a 180 deg ambi- 
guity. More advanced models that include up/downwind 
asymmetry due to effects such as Bragg scattering are 
currently under development. 

To estimate wind vectors from GPS observations we first 
determine the surface slope statistics in the model wrhich 
best fit the data. The second step is to use an ocean model 
to retrieve a wind speed and direction estimate from the 
slopes. Tw t o satellites are necessary to recover the direc- 
tional information lost in measuring the return from the 
entire annulus of the glistening zone with the DMR. In 
the current implementation, w ? e try to use satellites wuth 
relatively high elevation angles to insure that the spec- 
ular points are nearly co-located and that the elevation 



angle-dependent slope statistics are nearly the same for 
both satellites. A degenerate case occurs when the an- 
gle between the incident planes of the two GPS satellites 
is very close to 90 deg. In this situation only the mag- 
nitude of the angle between the incident plane and the 
up/downwind direction can be determined, and an addi- 
tional ambiguity exists in the wind direction. 

The correct wind speed is chosen by interpolating the es- 
timated slope statistics against the slope statistics from 
the Elfouhaily model [7]. Underlying the wave spectra 
and Z-V models are assumptions about the ocean condi- 
tions that affect the performance expected. In particular, 
the current model assumes well developed seas, with no 
swell. Coastal effects are not included. These simpli- 
fying assumptions limit the applicability of the current 
implementation under coastal and highly variable ocean 
conditions. 

INITIAL RESULTS 

On February 28 and March 11, 2000, Dr. Katzberg con- 
ducted two flight tests in a B200 aircraft off the Vir- 
gina coast to gather GPS reflection data, coincident with 
QuikSCAT overflights. Fig. 1 and Fig. 3 show the 
QuikSCAT wind fields and aircraft flight tracks for March 
11 and February 28, respectively. The February 28 map 
also shows the location of three buoys for which wind 
speed and direction data was logged by the National Data 
Buoy Center (NDBC). (Unfortunately, similar buoy data 
for March 11, 2000 has not been archived.) DMR data 
logged at 10 Hz were grouped into one-minute segments 
numbered sequentially along the flight tracks. Results of 
the GPS-based wind vector retrievals for each one minute 
sample on the two days are shown in Figs. 2 and 4. 

The March 11 QuikSCAT wind field shows very stable 
winds in the range from 8.5-10 m/s, due north (within +/- 
10 deg) over the entire aircraft flight path. GPS-derived 
wind direction estimates, shown in Fig. 2b are also con- 
sistently within about 5 deg of North. GPS-based wind 
speeds for the first 10 samples are quite low, reporting val- 
ues as low as 5 m/s as compared to the consistent 10 m/s 
QuikSCAT result. Beginning at sample 18, the GPS wind 
speed estimates move toward the QuikSCAT values with 
most of the points falling within the 2 m/s error range. 

This data set presents a very favorable situation for GPS 
versus QuikSCAT comparisons. The wind field is quite 
stable and does not originate over the land; thus the wave 
spectra model assumptions should be valid under these 
conditions and good agreement is to be expected. The 
discrepancy in wand speed may be related to a relatively 
large difference between the elevations of the tw'o satellites 
at the beginning of the flight. Initially one satellite is at 
an elevation of 44 deg while the second is at 63 deg. By 
the middle of the flight, the first has risen to 53 deg. We 
expect that a more accurate wund speed estimate can be 
obtained using only one of the two satellites or by modi- 
fying our algorithm to accommodate different elevations. 

The February 28 QuikSCAT wind field presented in Fig. 3 
shows quite a different situation with relatively low^ wind 


speeds (< 6 m/s) throughout, and wind directions that 
vary over 70 deg. Furthermore, winds along the flight 
path originate over the land areas. Surprisingly, GPS 
wind speed estimates fall within 2 m/s of the QuikSCAT 
estimates, despite the low wind speeds and a large differ- 
ence in satellite elevations. The GPS-derived wind direc- 
tion is very stable, in the range of 80-90 deg throughout 
the flight, not at all reflecting the rapid variations seen 
in the QuikSCAT estimates. There are number of possi- 
ble explanations for this, most important of which is the 
low wind speed, which may lead to unreliable direction 
estimates for both methods, A second possibility is the 
longer averaging time inherent in the GPS L-Band obser- 
vations and data grouping. Hourly data from all three 
buoys show T s consistent wind direction in the range of 90- 
110 deg over a two hour period surrounding the flight 
tests. A third consideration is related to the GPS satellite 
geometry. Throughout this flight, the azimuthal separa- 
tion of the tw r o GPS satellites w r as very close to 90 deg. 
As described previously, this degenerate case leads to an 
ambiguity in the wind direction related to the angle be- 
tween the w r ind and the incident plane. For a true wind 
direction of 90 deg, at the start of the experiment (Sample 
1), the ambiguity is such that the possible wind directions 
arc only 90 deg and 270 deg; however, for sample 30, the 
two satellite azimuths have shifted by 30 deg, resulting in 
possible wind directions of 90, 150, 270, and 330 deg. 

CONCLUSIONS 

The results presented here represent the first compar- 
isons between GPS-derived ocean surface wind vectors 
and QuikSCAT wind fields. We plan to continue these 
comparisons with additional flight experiments and anal- 
ysis under both favorable and challenging conditions with 
the goal of exploring the possibility of using GPS as a com- 
plimentary sensor to airborne and spaceborne scatterom- 
eters. The use of GPS for w r ind retrievals, is an emerging 
technology with considerable opportunity for further de- 
velopment of the theory and instrumentation, and signif- 
icant potential for contributing to our understanding of 
regional and global ocean surface conditions. 
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Figure 1: Quiver plot of QuikSCAT data for 11 March, 
2000. Also shown is the aircraft flight path. 
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Figure 2: Wind speed (a) and wind direction (b) estimates 
from GPS observations for 11 March, 2000 GPS data. 



Figure 3: Quiver plot of QuikSCAT data for 28 February, 
2000. Also shown are the aircraft flight path, and the 
locations of the Delaware Bay (DB), Chesapeake Light 
(CL), and Virginia Beach (VB) buoys. 
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Figure 4: Wind speed (a) and wind direction (b) estimates 
from GPS observations for 28 February, 2000 GPS data. 
Also shown on the plots are the Delaware Bay (DB) and 
Virginia Beach (VB) buoy estimates at the positions in 
the flight track closest to the buoy positions. 






